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ABSTRACT 

We study the radial migration of dust particles in accreting protostellar disks analogous to 
the primordial solar nebula. Our main objective is to determine the retention efficiency of dust 
particles which are the building blocks of the much larger planetesimals. This study takes account 
of the two dimensional (radial and normal) structure of the disk gas, including the effects of the 
variation in the gas velocity as a function of distance from the midplane. It is shown that the 
dust component of disks accretes slower than the gas component. At high altitude from the disk 
midplane (higher than a few disk scale heights), the gas rotates faster than particles because of 
the inward pressure gradient force, and its drag force causes particles to move outward in the 
radial direction. Viscous torque induces the gas within a scale height from the disk midplane 
to flow outward, carrying small (size < 100 /J,m at 10 AU) particles with it. Only particles at 
intermediate altitude or with sufficiently large sizes (> 1 mm at 10 AU) move inward. When 
the particles' radial velocities are averaged over the entire vertical direction, particles have a net 
inward flux. The magnitude of their radial motion depends on the particles' distance from the 
central star. At large distances, particles migrate inward with a velocity much faster than the gas 
accretion velocity. However, their inward velocity is reduced below that of the gas in the inner 
regions of the disk. The rate of velocity decrease is a function of the particles' size. While larger 
particles retain fast accretion velocity until they approach closer to the star, 10 /im particles 
have slower velocity than the gas in the most part of the disk (r < 100 AU). This differential 
migration of particles causes the size fractionation. Dust disks composed mostly of small particles 
(size < 10 /Ltm) accrete slower than gas disks, resulting in the increase in the dust-gas ratio during 
the gas accretion phase. If the gas disk has a steep radial density gradient or if dust particles 
sediment effectively to the disk midplane, the net vertically averaged flux of particles can be 
outward. In this case, the accretion of the dust component is prevented, leading to the formation 
of residual dust disks after their gas component is severely depleted. 

Subject headings: accretion, accretion disks — planetary systems: formation — solar system: formation 



1. Introduction 

Planets form in circumstellar disks. In the 
standard scenario, formation of earth-size planets 
or planetary cores occurs through coagulation of 
small dust particles (e.g., Weidenschilling & Cuzzi 
1993). Thus, the total amount, radial density dis- 
tribution, and size distribution of dust particles in 
the disks are the important initial conditions of 
planet formation. In many of the previous studies 
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of planet formation, the dust-gas ratio of circum- 
stellar disks is assumed to be similar to the solar 
value 10~^), and dust particles are considered 
to be well mixed with gas, i.e., the dust-gas ratio is 
constant throughout the disk (Hayashi, Nakazawa, 
& Nakagawa 1985). 

However, the motion of dust particles is differ- 
ent from that of gas. Dust particles have a radial 
motion which is induced by the gas drag force. 
Adachi, Hayashi, & Nakazawa (1976) and Weiden- 
schilling (1977) studied the motion of particles in 
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gas disks. Due to radially outward pressure gra- 
dient force, the rotation velocity of the gas is gen- 
erally slower than particles on nearly Keplerian 
circular orbits. Consequently, the gas drag force 
takes angular momentum from particles, resulting 
in their inward migration (Whipple 1972). At a 
few AU from the central star, the orbital decay 
time is estimated to be ~ 10^ yr for 100 /xm par- 
ticles and 10* yr for 1 cm particles. These times 
arc much shorter than the life time of gas disks 
(~ 10^ yr). Thus, dust component may evolve 
much faster than the gas in protostellar disks, in 
which case, the d\ist-gas ratio would decline. It 
appears that the initial conditions of planet for- 
mation need not to be protostellar disks in which 
dust particles are well mixed with the gas and the 
dust-gas ratio is constant. 

In this paper, we start a series of studies to 
determine the initial distribution of dust parti- 
cles in the context of planet formation. Here, we 
study the radial migration of dust particles in gas 
disks. We focus our attention on particles smaller 
than ~ 1 cm in order to consider evolution of dust 
disks before the initiation of planetesimal forma- 
tion. The migration of these small particles estab- 
lishes the initial density distribution of the dust 
component and sets the stage for planetesimal for- 
mation. In this phase of disk evolution, the gas 
component is considered to be turbulent and is ac- 
creting onto the star. The discussion in the above 
paragraph is based on the study of particle mi- 
gration under the assumption that particles have 
completely sedimented on the midplane of laminar 
disks. However, in turbulent disks, the sedimenta- 
tion of particles is prevented, because particles are 
stirred up to high altitude from the midplane by 
turbulent gas motion. In such disks, particles are 
distributed in the vertical direction, and the radial 
motion of particles depends on the distance from 
the midplane. Thus, there is a need for study- 
ing the radial migration of particles which reside 
above the midplane. 

There are two important factors causing the 
vertical variation in the migration velocity of par- 
ticles. The first is the variation of rotation velocity 
of the gas in the vertical direction. As mentioned 
above, the gas rotation differs from the Keplerian 
rotation because of the gas pressure gradient. As 
the gas density decreases with the distance from 
the midplane, the radial pressure gradient varies 



and under some circumstances even changes its 
sign. While at the midplane the pressure gradient 
force is outward, it is inward near the disk surface 
becaiisc the disk thickness increases with the ra- 
dius. The gas drag force on particles also varies 
with the height, resulting in a variation of particle 
migration velocities. In §3.1, we will see that par- 
ticles at high altitude, where the pressure gradient 
force is inward, flow outward. 

The second factor is the variation in the radial 
gas flow. Small particles (< 100 /im) are well cou- 
pled to the gas, so that they migrate with the gas 
flow, as discussed in §3.1 below. That is if the gas 
is accreted onto the star, the small particles would 
also be accreted, while if the gas flows outward, 
particle flow would also be outward. In accreting 
protostellar disks, the gas does not always flow ra- 
dially inward. Radial velocity of the gas varies 
with the distance from the midplane. Viscous 
stress can cause gas outflow near the midplane 
(Urpin 1984; Kley & Lin 1992). Rozyczka, Bo- 
denheimer, & Bell (1994) showed that the outflow 
occurs if the radial pressure gradient at the mid- 
plane is steep enough. The density (or pressure) 
gradient at the midplane is steeper than the av- 
erage spatial (or surface) density gradient, e.g., if 
the surface density varies as Eg oc and the disk 
thickness varies as hg oc r, the midplane density 
has a steeper variation corresponding to pg oc r~^. 
Viscous diffusion at the midplane causes the out- 
flow of the gas to reduce the steep density gra- 
dient. On the other hand, at high altitude from 
the midplane, the density gradient is shallow and 
the viscous torque causes the usual inward flow. 
Thus, in accretion disks, the outward flow at the 
midplane is sandwiched by the inflow at the disk 
surfaces (sec §2.2 for details). Because small parti- 
cles are carried by the gas flow, they are expected 
to flow outward near the midplane, and inward at 
higher altitude. 

Thus, if most of dust particles concentrate in 
the region where they migrate outward, i.e., near 
the disk surface or around the midplane, the dust 
component of the disk would accumulate in total 
mass and expand in size. In reality, the vertical 
spatial distribution of particles is determined by 
the degree of sedimentation, and their size distri- 
bution is regulated by the coagulation, condensa- 
tion, and sublimation processes. In this paper, we 
calculate the vertical distribution of particles in or- 
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der to identify the mass flux of the dust particles 
as functions of their size and their distance from 
the midplane and their host stars. The results of 
the present calculation will be used to study the 
long term dynamical evolution of the dust compo- 
nent of the disk. For the present task, we do not 
consider the effects of gas depiction, particles' size 
evolution or their feedback influence on the flow 
velocity of the gas. These effects will be consid- 
ered in a future investigation. 

There are several works on the particle mi- 
gration in turbulent disks. Stepinski & Valageas 
(1997) derived the vertically averaged velocity of 
migration for two limiting cases. One model 
is constructed for small particles which are well 
mixed with the gas and the other is for large parti- 
cles which are totally sedimented to the midplane. 
The velocity for intermediate sized particles is cal- 
culated by the interpolation of these limiting cases. 
Supulver & Lin (2000) investigated the evolution 
of particle orbits in turbulent disks numerically. 
This paper gives analytical expression of the parti- 
cle velocity for any size of particles, which are esti- 
mated from interpolation by Stepinski & Valageas 
(1997). Several studies showed that vortices or 
turbulent eddies can trap particles within some 
size range (Barge & Sommeria 1995; Cuzzi, Dobro- 
volskis, & Hogan 1996; Tanga et al. 1996; Klahr 
& Henning 1997; Cuzzi et al. 2001, see however 
Hodgson & Brandenburg 1998). Particle trapping 
by eddies may be important for local collection 
of size-sorted particles. In this paper, turbulent 
eddies are treated just as a source of the particle 
diffusion and of the gas viscosity. The effect of par- 
ticle trapping on their coagulation is not discussed 
here. 

The plan of the paper is as follows. In §2, the 
vertical variation of gas flow is derived. In §3, we 
describe how the radial velocity of dust particles 
varies in the vertical direction, and then calculate 
the net radial migration velocity. In §4, we dis- 
cuss the steady distribution of particles assuming 
no particle growth. We show that the size frac- 
tionation of particles occurs as a result of radial 
migration, and that the distribution of particles 
differs from that of the gas. 



2. Viscous Flows of Gas Disks 
2.1. Rotation Law of the Gas 

Gas disks rotate with a slightly different veloc- 
ity from the Keplerian velocity. The rotation ve- 
locity is determined by the balance between the 
stellar gravity, the centrifugal force, and the gas 
pressure gradient. In cylindrical coordinates (r. z), 
the balance of the forces described by the momen- 
tum equation in the r-direction is given by 
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where Qg is the angular rotation velocity of the 
gas, G is the gravitational constant, M is the stel- 
lar mass, pg is the gas density, and Pg is the gas 
pressure. The density distribution in the vertical 
direction is determined by the balance between the 
^-components of the stellar gravity and the pres- 
sure gradient, which is written as 
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We assume that the disk is isothermal in the ver- 
tical direction. Neglecting (zjr)^ and higher order 
terms, the vertical density distribution is obtained 
as 

(3) 



Pg(r, z) = pg{r, 0) exp (^-^ 
The disk scale height hg is given by 
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^^K.mid 

where c is the isothermal sound speed, and 
^^K,mid = y/GM/r^ is the Keplerian angular ve- 
locity at the midplane. The radial variations of 
physical quantities are assumed to have power law 
forms, such that 

Pg{r,z) = porluexp (^-^^ , 



2^9 



O'AU ' 

„(9+3)/2 



(5) 



where the subscript "0" means the values at 1 AU, 
tau is the radius in the unit of AU, and power law 
indices p and q are usually negative. The surface 
density is 
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where ps = P + (<Z + 3)/2. To derive the rotation 
law of the gas, we retain {z/r)"^ and lower order 
terms in equation (1), then obtain 



i^g{r,z) = fiK.mid 
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The difference of the gas rotation from the Kepler 
rotation is order of {hg/rY, and the gas at high 
altitude rotates slower than the gas at the mid- 
plane. 

2.2. Radial Velocity of the Gas Flow 

In viscous disks, angular momentum is trans- 
fered from the inner part to the outer part of 
the disk under the action of viscous stress. The 
azimuthal component of the momentum equation 
implies 
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where Vr,g and v^^g are the r- and z-components of 
the gas velocity, respectively, and v is the kinetic 
viscosity. The right hand side of equation (8) rep- 
resents the torque exerted on an annulus with unit 
width in the r- and 2r-dircctions. The left hand side 
is the variation of the angular momentum of the 
annulus accompanying with the motion. We as- 
sume that the angular velocity is always adjusted 
to the equilibrium (eq. [7]) on a dynamical time 
scale. Since molcciilar viscosity is too ineffective, 
it is customary to assume that angular momentum 
is primarily transported by the Reynolds stress in- 
duced by the turbulent motion of the gas. Follow- 
ing conventional practice, we adopt an ad hoc but 
simple to use a prescription for turbulent viscosity 
(Shakura & Sunyaev 1973), such that 



V = achg = acohor'j^^^ . 
The mass conservation of the gas is 
dp„ Id, , d , , 
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In accretion disks considered here, the second term 
is the same order of the third term. Then, we 
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Thus, Vz.gd{r'^Vlg) / dz is factor {hg/rY smaller 
than Vr,gd{r^ng)/dr and can be neglected in equa- 
tion (8). Using power law expressions (5) and (9) 
and retaining terms of order (hg/rY, the radial 
velocity is reduced to 
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In the following discussions, we adopt these val- 
ues as being representative of a typical model: 
M = IMq, po = 2.83 X 10-1° g cm'^, ho = 
3.33 X 10-2 AU, p = -2.25, q = -0.5, and 
a = 10"^. These values correspond to a min- 
imum mass solar nebula whose surface density 
distribution has a power law form with index 
Ps = P + {q + 3)/2 = —1.0 and has a gas mass 
2.5 X IO-^Mq inside lOOAU. 

Figure 1 shows the radial velocity of the gas. 
The radial velocity is outward (positive) near the 
midplane, while it is inward (negative) above the 
height 0.73hg. In the standard model (q = —0.5), 
the radial velocity is independent of the distance 
from the star, as seen from equation (11). At the 
midplane, the radial gradient of the gas density 
is so steep {p = —2.25) that the gas receives more 
torque from the inner disk than the torque it exerts 
on the outer disk. This radial dependence means 
that viscous stress acts to move the gas outward 
in such a steep density gradient. The condition for 
the outward gas flow at the midplane is 



p+-q<-2 



(12) 



which is satisfied in many accretion disk models. 
The rotation of the gas is slower at higher altitude 
as seen in equation (7). Thus, the gas at the mid- 
plane loses its angular momentum through the ac- 
tion of viscous stress it has with the gas at higher 
altitude. However, in the standard model, this 
angular momentum loss is smaller than the an- 
gular momentum gain through the viscous stress 
in the r-dircction. The outflow zone around the 
midplane shrinks as the radial density gradient is 
reduced (p — ^ 0), and when the condition (12) is 
violated, gas at all altitudes flows inward. The 
radial density gradient is shallower at higher alti- 
tude, because the disk scale height increases with 
the radius. Thus, at high altitude, the gas loses 
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more angular momentum through exerting torque 
on the outer disk than it is receiving from the inner 
disk. For such gas, the presence of viscous torque 
induces it to flow towards the star. 

The net radial velocity of the gas is given by 
averaging in the vertical direction: 



^ I Vr,gPgdz 
J-oo 
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The net accretion of the gas toward the star occurs 
if 

Ps + q^P+^{q + l)> -2 . (14) 

We may consider two models in which a steady 
state is achieved. One model is the case in which 
the net radial velocity is zero, i.e., Ps+q = —2. For 
example, if Ps = —1.5 and q = —0.5, there is no 
net accretion of the disk. A more realistic model 
is the case in which the mass flux of the gas is con- 
stant with the radius, i.e., rT,g{vr,g) oc 
is constant. We take this as our standard model, 
i.e., ps = —1.0 and q = —0.5. This surface den- 
sity profile is expected to arise after the disk has 
undergone a period of initial viscous evolution. 

3. Radial Flov^f of Dust Particles 

3.1. Vertical Dependence of Radial Veloc- 
ity 

The azimuthal velocities of dust particles are 
different from that of the gas. The resulting gas 
drag force transfers angular momentum between 
the particles and the gas and moves particles in 
the radial direction. If there is no gas drag force, 
particles would orbit with the Keplerian angular 
velocity, which is approximated as 
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(15) 



We express the deviation of the angular velocity 
of the gas from the Keplerian angular velocity as 



fig = fiK(l - V) 



(16) 



Prom equation (1) it is seen that 

T] = —{rfl^pg)~^ dPg /dr is the ratio of the gas 



pressure gradient to the stellar gravity in the radial 
direction. From equations (7), (15), and (16), rj is 
written within the order of (hg/r)'^ as. 
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Note that the pressure gradient force is outward {r] 
is positive) around the midplane, while it is inward 
[r] is negative) where \z\ > ^/—2{p + q)/{q + •i)hg. 
In the standard model, ry changes sign at ^ « 
1.5hg. Dust particles near the midplane rotate 
faster than the gas, and particles at high altitude 
rotate slower than the gas. 

The equations of motion of a particle are 



dVr,d _ '^e,d --,2 fiK,inid 
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{Vr,d - Vr,g) , (18) 



(19) 



where Vr and are the r- and 9- components of 
the velocity, respectively, with the subscripts "(7" 
and "d" distinguishing gas and dust, and VK.mid = 
f'^K,mid is the Keplerian velocity at the midplane. 
The gas drag force is expressed through the non- 
dimensional stopping time Tg (normalized by the 
Kepler time at the midplane). We do not solve the 
equation of motion in the z-direction. Instead, we 
simply assume Vz,d = 0. When an equilibrium in 
the vertical dust distribution is achieved, the dust 
sedimentation to the midplane is balanced by the 
diffusion of particles, as discussed below in §3.2. 
The vertical velocity of a particle is zero when time 
averaged. 

The mean free path of gas molecules is larger 
than 1 cm for r > 1 AU in our models (Naka- 
gawa, Sekiya, & Hayashi 1986). In this paper, we 
consider particles smaller than 1 cm and use Ep- 
stein's gas drag law. Then, the non-dimensional 
stopping time T^, as given by Takeuchi & Arty- 
mowicz (2001), is 



T, = 



PpSt^K.mid 



(20) 



PgTVT 

where pp is the particle internal density, s is the 
particle radius, and the mean thermal velocity is 



We take the particle internal den- 



tin Takeuchi & Artymowicz (2001), v^p is defined to be 4/3 
times the mean thermal velocity. This definition causes a 
factor 4/3 difl^erence between equation (20) and their equa- 
tion (10). 
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sity as Pp = 1.25 g cm ^. 

For particles smaller than 1 cm, the non- 
dimensional stopping time is much smaller than 
unity through most of the disk. Figure 2 shows 
the radii of particles whose stopping time is unity. 
At the midplane (see the solid line), the non- 
dimensional stopping time is smaller than unity 
for particles of size s < 1 cm at r < 100 AU. At 
z = 2hg (the dashed line), although the gas den- 
sity is lower than at the midplane and the stopping 
time is longer, particles of size s <1 mm still have 
the non-dimensional stopping time smaller than 
unity. These particles are well coupled to the gas 
and their angular velocity is similar to the gas 
angular velocity (7). Only particles with s > 1 cm 
which are located at high altitude and r > 100 AU 
become decoupled from the gas. 

We assume that motions of both the gas 
and the particles are close to Keplerian, i.e., 

ve^g ~ ve,d ~ "K,mid, and that d{rvo,d)/dt » 
Vr,dd{rvK,mid) / dr = Wr,dfK,mid/2. Then, from 
equation (19), we have 

ve,d - ve,g = -^TsVr,d ■ (21) 

Using equation (16) and neglecting terms of order 
{hg/r)'^ and higher, equation (18) is reduced to 
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{Vr,d - Vr,g) . (22) 



The left hand side is order of v^Jr and is ne- 
glected if Vr^d ^ c. Substituting equation (21) 
into equation (22), we find the radial velocity of 
the particle to be 



??l^K,mid 



(23) 



In the above derivation of the particle's radial ve- 
locity, we assume that in the z-dircction the par- 
ticle sedimentation is balanced by the turbulent 
diffusion, but that in the r-direction the turbulent 
effects are neglected. 

For particles well coupled with the gas (T^ <C 
1), the radial velocity reduces to 



vr,g 



■ Vr,drift 



(24) 



where Vr^AnMt — ^'yTsWK.mid is the relative velocity 
from the gas. Using equations (5), (17), and (20), 



I'r.drift = 27r 



AU 
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X rX+'/'r,,wexp AU yr'^ , (25) 

where Tsjinid is the non-dimensional stopping time 
at the midplane. The radial drift velocity increases 
exponentially with the height. Because the parti- 
cles are strongly coupled with the gas, the gas drag 
force suppresses their drift velocity. At the mid- 
plane where the gas density is highest, the sup- 
pression of the drift is most effective. The drift 
velocity at the midplane is 



^r.drift.mid 



[t^ hoS Pp. , ■, -p+q/2-l .j-r ^1 



(26) 

which agrees with the derivations of Adachi et al. 
(1976) and WeidenschiUing (1977). 

For large particles decoupled from the gas 
{Tg ^ 1), the radial motion of the gas does 
not affect the particles' velocity, so that Vr,d = 
— ?7T~^t;K,mid- In this paper, we do not consider 
such large particles. 

Figure 3 shows the radial velocity Vr^d of parti- 
cles of s = 10 /um, 100 /im, and 1 mm at 10 AU 
in the standard model disk. Near the midplane, 
the drift velocity Wdrift (relative to the gas veloc- 
ity shown by the dotted line) of 10 /im particles 
is small. The particles move outward almost to- 
gether with the gas. As the altitude increases 
and the gas velocity decreases, the radial veloc- 
ity of the particles also decreases, then becomes 
negative. At 2; ~ 1.5/ig where r] changes its sign, 
the radial pressure gradient vanishes and the gas 
rotates with the Keplerian velocity. The parti- 
cles co-rotate with the gas at that location, and 
have the same radial velocity as the gas. The 
drift velocity changes its sign to become positive. 
At high altitude {z > 2hg) where the gas drag 
is weak, the particles' drift velocity begins to in- 
crease rapidly, and then the radial velocity be- 
comes positive again at z > 3hg. The behavior 
of the radial velocity of 100 /im particles is quali- 
tatively similar to that of 10 p,m. particles, i.e., it is 
outward at the midplane, inward at intermediate 
altitude, and outward again at high altitude. The 
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coupling of 100 fim particles to the gas is weaker, 
so that their radial velocity is much different from 
the gas velocity even at the midplane. However, 
the drift velocity at the midplane is still slightly 
smaller than the gas outflow velocity and the par- 
ticles move outward. Particles of 1 mm move in- 
ward even at the midplane, because the inward 
drift velocity of such large particles is larger than 
the outflow velocity of the gas. These large parti- 
cles also move outward at high altitude where the 
gas rotates faster than the particles. 

3.2. Vertical Distribution of Particles 

At the beginning stage of star formation from 
molecular cloud cores, dust particles in circum- 
stellar disks may be well mixed with the gas and 
the dust-gas ratio may be uniform throughout the 
disk. After the termination of gas infall onto the 
disks, the disks become to be hydrostatic and par- 
ticles begin to sediment toward the midplane be- 
cause of the z-component of the stellar gravity. 
The particles reach the terminal velocity, where 
the gravity and the gas drag balance with each 
other, 

Vz,d = -i^KMdTsZ . (27) 

The time scale of the sedimentation is tsed ~ 
z/vz.d ^ T'r^^K^- If this time scale is much 
smaller than the radial migration time scale, tr ~ 
r/vr,d ~ (a + T,)-^ihg/r)-^n-\ i.e., if T,/a > 
(hg/rY, then the particles sediment before the 
large migration in the r-direction. (At high alti- 
tudes or around the midplane, particles may drift 
to the opposite direction to the gas, as seen in 
Fig. 3. Some particles with ~ a may have the 
radial migration time scale tr much larger than 
the above estimate. Even for such particles, the 
condition Tg/a 3> {hg/r)^ for the fast sedimenta- 
tion is still valid.) For example, at 10 AU in a 
disk with a = 10~^ and (hg/r)"^ ~ 10"^, parti- 
cles larger than 0.2 /xm sediment to the midplane 
without large radial movement. 

We consider disks at later stages, such as the 
T Tauri stage. The disks are still turbulent, and 
the turbulent motion of the gas stirs dust particles 
up to high altitude to prevent dust sedimentation. 
An equilibrium distribution of particles in the ver- 
tical direction is achieved by the balance between 
the sedimentation and the diffusion due to the tur- 
bulent gas. 



The turbulent diffusion is modeled in an anal- 
ogy of molecular diffusion, provided that the par- 
ticles are considered to be the passive tracers of 
fluid, i.e., if the particles have no influence on the 
gas motion and have the similar velocity of the 
surrounding gas [see e.g., chapter 10 in Monin & 
Yaglom (1971) and section 3.5.1 in MorflU (1985)]. 
The equation of continuity is written as 

g-^Pd + V-{pdVd + j) = 0, (28) 

where pd is the particle density, and f d is the parti- 
cle velocity. The diffusive mass flux j is estimated 

where the Schmidt number Sc represents the 
strength of coupling between the particles and 
the gas. For small particles, Sc approaches unity 
and the particles have the same diffusivity as the 
gas, while it becomes infinite for large particles. 
For intermediate particle sizes, Sc can be as small 
as 0.1, which means that particle diffusion occurs 
effectively (see e.g.. Fig. 1 in Cuzzi et al. 1993). 
In our standard model, we use Sc = 1. If Sc is 
not unity, the particles are not the passive tracers 
which completely follow the fluid motion. In this 
case, the formulation according to the analogy 
of molecular diffusion may be inappropriate. In 
addition, if the velocity of sedimentation, Vz^d, is 
comparable to or larger than the turbulent veloc- 
ity, the particle could not be the passive tracers. 
The estimate of the diffusive mass flxix by equa- 
tion (29) should be considered just as a "gradient 
diffusion hypothesis." 

In a steady state of axisymmetrical disks, d/dt 
and d/d0 in equation (28) is zero. In addition, for 
particles satisfying the condition Tg/a ^ (hg/r)^ 
(i.e., for particles sedimenting fast without large 
radial migration), we see that d{pdVz^d) / dz ^ 
d{pdVr,d)/dr. Therefore, the mass flux in the z- 
direction must be zero, i.e., 

P..M-P4^(g)=0. (30) 

This equation is the same as the one derived by 
Dubrulle, MorflU, & Sterzik (1995). Solving equa- 
tion (30) with equations (3), (4), (9), (20), and 



7 



(27) gives the particle density 



Pd{r,z) = pd{r,0)ex.p 



ScT, 



5, mid 



2hl 



a 



exp 



2hl 



(31) 

The first term in the exponential comes from the 

gas distribution and the second term represents 
the sedimentation. In the limit of tight coupling of 
the particles and the gas (Tg^mid = 0), the particle 
distribution is the same as that of the gas. The 
surface mass density of dust particles is 



Mr) 



+ 00 



Pd{r, z)dz 



(32) 



Figure 4a shows variation of the dust-gas ratio 
in the vertical direction at 10 AU in the standard 
model disk. We see the sedimentation of particles. 
^ Large particles (> 1 mm) sediment around the 
midplane, while small particles (< 10 /xm) spread 
over a few scale height of the gas disk. These small 
particles are stirred up to high altitude by the tur- 
bulence of the gas. The gas drag force becomes 
weaker at higher altitude where the gas density 
decreases, so the gas cannot sustain the dust parti- 
cles there. Thus, the particle density drops rapidly 
at some altitude and the dust disk obtains a rela- 
tively sharp surface (e.g., at z ~ 2.5/ig for 10 /zm 
particles). The sedimentation of particles is more 
effective in the outer part of the disk, farther away 
from the star (see Fig. 46). For example, 100 /zm 
particles at 1 AU spread over the gas disk, while 
the particles at 100 AU concentrate aroimd the 
midplane. At the outer part of the disk, the gas 
density is lower, so that the gas drag is weaker 
and turbulence cannot loft the particles to high 
altitudes. 

Dust particles are well mixed with the gas if 



ScT, ,mid I Z 

I exp 777"^ — 1 



a 



2hl 



i.e., below the height 



z 

hn 



2 log 



ScT, 



+ 1 



5, mid 



1/2 



(33) 



(34) 



^We obtained slightly thinner distribution of particles than 

that by DubruUc ct al. (1995) which is shown in their Fig- 
ure 3. They assumed that the stopping time of particles, 
which is inversely proportional to the gas density, to be con- 
stant and used the value at the midplane. This assumption 
causes the dust disk to puff out. 



In the standard model, particles smaller than 

s «; 20 r^u Mm (35) 

uniformly distribute with regard to the gas under 
the height z = 3hg. 

3.3. Net Radial Velocity 

As discussed in §3.1, the radial velocity of par- 
ticles varies with the altitude. In this subsection, 
we discuss the net radial velocity averaged in the 
vertical direction. 

3.3.1. Particles Well Mixed with the Gas 

At the beginning stage of formation of circum- 
stellar disks, small particles would be distributed 
uniformly in the disk gas. Before the particles have 
undergone enough sedimentation, their net radial 
velocity averaged in the vertical direction is 



1 f+°° 

J-oo 



pgVr,ddz . 



(36) 



Figure 5 shows the net radial velocities of small 
particles. The velocity of particles approaches 
the gas velocity as the distance from the star 
becomes smaller. Near the star, the gas den- 
sity is high enough to induce almost all particles 
to move together with the gas. As the distance 
from the star increases, the gas density decreases 
and particles at high altitude begin to drift out- 
ward from the gas. At sufficiently large distance 
from the star, the net radial velocity can be pos- 
itive, i.e., particles move outward. However, be- 
fore they move over a large distance in the ra- 
dial direction, they also sediment to the midplane. 
The time scale of the outward radial motion is 
tout ~ r/ {vr,d) ~ l/{r}Ts^K) (for outflowing parti- 
cles, \vr,g\ is smaller than |wr,drift| and is neglected 
in the estimate of tout); while the sedimentation 
time scale is tg^d ~ z/v^^d ~ 1/(Is^k)- Because 
?7 <C 1, tout ^ ^sod- The outward motion of those 
particles which are well mixed with the gas may 
proceed shortly after the formation of circumstel- 
lar disks, but it does not significantly modify the 
radial distribution of dust particles. The evolution 
of radial distribution occurs primarily through the 
motion of sedimented particles. 
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3.3.2. Removal of Outflow at High Altitude by 

the Sedimentation 

The net radial velocity of sedimented particles 

is 

1 /■+°° 

{Vr,d) = PdVr,ddz . (37) 

J-oo 

The function PdVr.dl^d-, which is proportional to 
the mass flux, is shown in Figure 6. Because of the 
sedimentation to the midplanc, which is efficient 
for larger particles, the mass flux is dominated 
by the particles near the midplane. For example, 
the mass flux of 10 /im particles is mainly car- 
ried by particles moving outward around the mid- 
plane (|^;| < 0.7 kg) and by those moving inward 
at intermediate altitude (0.7hg < \z\ < 2.9hg). 
Although particles at high altitude {\z\ > 2.9/ig) 
move outward, their contribution to the mass flux 
is negligibly small. Very large particles, for exam- 
ple 1 mm particles, move inward even at the mid- 
plane. Thus, the sedimentation causes the vast 
majority of particles to move inward. 

Sedimentation effectively removes outflowing 
particles at high altitude. This flow pattern is seen 
as follows. The particles flowing outward have a 
drift velocity larger than the gas inflow velocity, 
i.e., |fr,drift| > \vr,g\- Prom equatious (11) and 
(25), it is seen that'a-^Ts^mid exp(2;V2/ip > 1 for 
such particles. However, from equation (31), we 
also see that the density of these particles is as 
small as pd{z) / pd{0) ^ e^^. Therefore, very few 
particles remain in the outflow region at high al- 
titude. 

If the sedimentation is so effective that most 
particles concentrate around the midplane (l^l < 
kg), e.g., for particles larger than 1 mm at 10 AU 
(see Fig. 4a), the inward drift velocity at the mid- 
plane is larger than the gas outflow velocity. Thus, 
such particles flow inward. Again from equation 
(31), we see that Ts^mid/ce > 1 for such particles. 
Then, it is seen from equations (11) and (25) that 
br,drift| ?J |'''r,g| at the midplane. This inequality 
implies that when particles grow up suSic;iently 
large sizes to sediment toward and to concentrate 
around the midplane, they become decoupled from 
the outflow motion of the gas. If the sedimenta- 
tion is so efflcient that all particles concentrate in a 
thin layer at the midplane (Tj^mid/a ^ 1), the ra- 
dial motion of the particles would not affected by 
the radial motion of the gas. In this limit, the par- 



ticles' motion is similar to that deduced by Adachi 
et al. (1976) and Weidenschilling (1977). 

3.3.3. Net Velocity of Sedimented Particles 

Figure 7a shows the net radial velocities of sed- 
imented particles. These net radial velocities are 
negative for particles of all sizes. Particles rapidly 
move inward when they are at large distances from 
the star. In the outer regions of the disk, par- 
ticles sediment and concentrate at the midplane, 
and their net radial velocity, which is much faster 
than the gas outflow velocity Vr^g, is approximated 
by the drift velocity at the midplane (eq. [26]). 
Their inward velocity decreases as they approach 
the star, because the suppression of inward veloc- 
ity by the gas drag becomes stronger. The time 
scale of their orbital decay is r/{vr^d) oc r'^+P~i/'^ . 
When they approach the location where the gas 
density is dense enough to make Tg^-^i^/a < 1, 
the inward drift velocity becomes smaller than the 
gas outflow velocity at the midplane. The parti- 
cles begin to move with the gas flow. At the same 
time, particles are spreaded over more than the 
disk scale height. The particles' net radial veloc- 
ity approaches the net gas velocity (the dashed line 
in Fig. 7a). However, because the particle distri- 
bution concentrates slightly more to the midplane 
than the gas distribution, the number of particles 
riding on the outflowing gas around the midplane 
is larger than in the case where the particle dis- 
tribution is same as the gas. The enhancement of 
outflowing particles retards their net inflow veloc- 
ity to values below that of the gas. At the inner 
part of the disk, the net inward velocity of particles 
is always slower than the gas velocity. Near the 
star, the particles mix more with the gas and their 
motion converges with that of the gas. The par- 
ticles' inward velocity has a minimum magnitude, 
{vr,d) = 5-2 X 10~^ AU/yr, which is about half of 
the gas velocity. Particles smaller than 10 /xm have 
a slower inward velocity than the gas throughout 
nearly the entire disk (r < 100 AU), while parti- 
cles larger than 1 mm have slower velocity only in 
the innermost part of the disk (r < 1 AU). The 
difference in the inward velocity induces the size 
fractionation of particles, as discussed in §4.1 be- 
low. The dust-gas ratio may increase during the 
viscous evolution of disks, because the accretion 
velocity of small particles is slower than that of 
the gas. 
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Small particles which satisfy the condition (35) 
are mixed well with the gas. The net radial ve- 
locities of such small particles may be calculated 
by equation (36), even after larger particles have 
undergone the sedimentation. It is seen from com- 
parison between Figures 5 and 7o that, e.g., for 
10 /im particles, the radial velocities arc about the 
same in both figures for r < 1 AU, while they de- 
viate significantly between these figures for r > 10 
AU. As shown in Figure 5, the radial velocity of 
10 /um particles at r > 10 AU is outward when 
the particles are mixed with the gas, i.e., in the 
very first stage of disk formation. However, af- 
ter the particle sedimentation, the radial velocity 
becomes inward even at r > 10 AU, as shown in 
Figure 7a. It is concluded that the particle sedi- 
mentation suppresses the outward flow of particles 
in the standard model. 

3.3.4. Various Models 

In Figure 7b, we show the case with ps = —0.5 
where the surface density profile of the gas is flat- 
ter than in the standard model. The velocity pro- 
files are qualitatively similar to those in the stan- 
dard model. Particles of 10 /xm move slower than 
the gas throughout nearly the entire disk, while 
1 mm particles rapidly migrate inward. 

In Figure 7c, the case with a steeper surface 
density gradient is shown. In this model, we adopt 
Ps = —1.3 which satisfies the condition (14) for in- 
ward gas accretion flow. In this case, we see that 
the net radial velocity of particles becomes posi- 
tive at some locations. In these regions, the accre- 
tion of particles is prevented. The particles flowing 
inward from large distances terminate their migra- 
tion at the location where the net radial velocity 
becomes zero, and then accumulate at that par- 
ticular orbital radius. Particles located just inside 
this critical radius flow outward and accumulate 
there also. Particles at the innermost part of the 
disk flow inward onto the star. Size fractiona- 
tion occurs because the location of the accumu- 
lation depends on the particle size. The accumu- 
lation continues as particles drift radially inward 
from large radii until particles in the outer disk 
are signiflcantly depleted or the number density 
of accumulated particles becomes so high that the 
removal of particles through coagulation or col- 
lisional destruction becomes eflicient. Although 
accretion of dust particles terminates at some lo- 



cations, the gas disk continues to accrete onto the 
star. The dust-gas ratio increases as the accretion 
of the gas proceeds. 

3.3.5. Self- similarity of the Velocity Profile 

The functions of net radial velocities shown in 
Figure 7 for various sizes have self-similar forms. 

From equations (11), (24), (25), (31), (32), and 
(37), the net radial velocity can be written as 



AU 



(38) 



where is a function of p, q, Sc, and Tg^^i^/a, 



F(p,q,Sc; 



Sjinid 

a 

r+00 



/_r{-3p-2g-6-(5g + 9). 



/21 T,, 



X exp 

+00 

exp 



■ exp z' 



-z'^ - (exp0'2 - 1)] dz' 



^,2 ScT,,„i. 



(expz'2-1) dz'^ 



and z' = z/V2hg. If we vary the properties of par- 
ticles, keeping Tg^mid/ct to be constant, the fimc- 
tion F yields the same value. For example, when 
we vary the particle size s s' , then transform 
the orbital radius as r ^ r' = {.s' /s)'^/^P+^i+^y'^h, 
we have the same value of rnid/a and the same 
functional form of F. In our models, q — —1/2, so 
{i'r,d) depends on r only through the function F. 
In this case, the profiles of the net radial velocities 
for different sizes are obtained by the transforma- 
tion in the r-direction from one profile. In Figure 
8, the function F is plotted for various values of 
the Schmidt number Sc. 

3.3.6. Various Schmidt Numbers 

In the standard model, the net radial veloc- 
ity of particles is inward everywhere in the disk. 
This flow pattern is partially due to the removal 
of particles, through sedimentation, from high al- 
titudes where they would flow outward. For par- 
ticles which are sufficiently large to concentrate 
near the midplane, the magnitude of their inward 
drift velocity is larger than that of the outflow- 
ing gas. Thus, the strength of the sedimentation 
is always moderate enough to remove outflowing 
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particles. The sedimentation is controlled by cou- 
pling of particles to the turbulent motion of the 
gas. If coupling is relatively strong (Sc is smaller 
than unity) , the sedimentation would be less effec- 
tive and particles would mix more completely to 
the gas, while if it is relatively weak (Sc is larger), 
particles would concentrate more to the midplane. 
Usually Sc is assumed to be unity, but the actual 
value may be different from unity. Some exper- 
imental data arc plotted in Figure 1 in Cuzzi et 
al. (1993) which show Sc as small as 0.1 when 
the Stokes number St is about 0.01 (the Stokes 
number is considered to be the same order of the 
non-dimensional stopping time Tg in our models). 
For large particles (St ~ Tg > 1), Sc increases with 
the particle size. 

Figure 7d shows the net radial velocities in the 
case where the sedimentation is efficient (Sc = 10). 
The net radial velocity becomes outward at its 
peak. Thus, an accumulation of particles and an 
increase in the dust-gas ratio would occur. Be- 
cause of the strong sedimentation, the dust disk is 
much thinner than in the standard model (Sc = 1) 
and concentrates more to the midplane where the 
gas flows outward. In the outer part of the disk 
(for example r > 10 AU for 100 /xm particles), the 
inward drift velocity of particles (relative to the 
gas) at the midplane is faster than the outward gas 
velocity. As the particles drift in, the inward veloc- 
ity decreases, and at some distances from the star 
(r ~ 10 AU for 100 /xm particles), the inward drift 
velocity at the midplane becomes smaller than the 
outward gas velocity, and the particles around the 
midplane flow outward. At such distances, the 
dust disk is still concentrated around the mid- 
plane, and the net radial velocity is dominated 
by the particles around the midplane and is out- 
ward. As the distance from the star decreases, the 
particles mix more with the gas and the concen- 
tration to the midplane becomes weaker. The net 
radial velocity approaches the gas velocity, which 
is inward, at the innermost part of the disk. 

If the Schmidt number Sc is less than unity, the 
sedimentation of particles is less efficient. There- 
fore, the particles are always well mixed with the 
gas at where the inward drift velocity of parti- 
cles is smaller than the gas outfiow velocity (i.e., 
at where the particle outflow at the midplane oc- 
curs). Thus, the net radial velocity of particles is 
inward everywhere. 



Figure 8 shows the function F for various values 
of the Schmidt number, Sc = 0.1, 1, and 10. The 
net radial velocity has the same sign as F (see eq. 
[38]). It is seen that the outflow of dust particles 
occurs if the sedimentation is efficient (Sc > 10). 

4. Discussion and Summary 

4.1. Steady Density Distribution of Dust 
Particles 

Because the inward velocity of particles de- 
pends on their size, particles of different sizes ac- 
cumulate at different locations. In this subsection, 
we discuss the density distribution of dust parti- 
cles as a consequence of their radial flow and show 
how their size fractionation may occur. We adopt 
the models in which the net velocity of particles 
is always inward (Sc = 1 and ps > —1.0). We as- 
sume that the distribution of particles approaches 
a steady state after a brief stage of initial evolu- 
tion. When a steady state is achieved, the mass 
flux of dust particles becomes constant in the ra- 
dial direction. We calculate the mass flux of par- 
ticles as 



dMd 

ds 



as 



(40) 



where Ma and are the mass flux and the sur- 
face density, respectively, of particles smaller than 
size s. For simplicity, we neglect three physical 
processes. First, the evolution of the particle size 
through the coagulation, collisional destruction, 
condensation, and sublimation is neglected, i.e., 
we assume the mass flux dM^/ds for each size 
range is constant in the radial direction. Second, 
in turbiflcnt disks, the particles' mass flux comes 
not only from the mean flow with an average veloc- 
ity {vr,d.) but also from the turbulent diffusion of 
particles, which appears as j in equation (28) for 
example. We neglect the mass flux from the tur- 
bulent diffusion. Third, we assume the structure 
of the gas component of the disk is entirely deter- 
mined by the gas itself. We neglect the feedback 
drag induced by the particles on the gas. In the 
limit that the spatial density of the dust compo- 
nent becomes comparable to that of the gas near 
the midplane, this effect would speed up the az- 
imuthal velocity of the gas to the Keplcrian value 
and quench the radial migration of the dust (Cuzzi 
et al. 1993). Such a dust concentration requires 
substantial sedimentation of relatively large parti- 
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clcs. Contributions from all of these factors will be 
investigated in future grain-evolution calculations. 
Here, we adopt the simplest assumptions to focus 
on showing the size fractionation of particles. The 
mass flux of particles in all size range is 

M,,au = £°''^ds, (41) 

where Smin and Smax are the minimum and max- 
imum sizes of particles, respectively. The surface 
density in the steady state is calculated from equa- 
tion (40) with given M^. Figure 9a shows the sur- 
face density of the dust component composed of 
single size particles. Particles of different sizes 
have different density profiles. If the dust com- 
ponent is composed of relatively large particles, it 
would be concentrate at the inner region of the 
disk. Thus, as particles grow in size, their sur- 
face density distribution becomes more centrally 
concentrated. Figure 9b shows the surface den- 
sity, assuming the size distribution of particles is 
a power law with index —3.5, Smin = 0.1 fim, and 
Smax = 10 cm. The surface density distribution 
of the dust particles is different from that of the 
gas. The power law index of the dust distribu- 
tion Tid is about —1.5 for a gas disk with index 
Ps = —1.0 (—1.2 for a gas disk with = —0.5). 
The implied power law index, —1.5, is similar to 
the value anticipated from the present mass distri- 
bution of planets in the solar system (Hayashi et 
al. 1985). However, note that the density distri- 
butions in Figure 9 are derived assuming no size 
evolution of particles, and no turbulent diffusion 
in the radial direction. The growth of particles 
during the radial flow adds a source term in the 
equation of continuity. The evolution of the dust 
density profile should be investigated further by 
taking particle growth and turbulent diffusion into 
account. 

4.2. Evolution of the Dust-gas Ratio 

The accretion velocity of dust particles is dif- 
ferent from that of the gas. This difference causes 
evolution of the dust-gas ratio. In the standard 
model, 10 /um particles between 10 and 100 AU 
have an inflow velocity which is about half of the 
gas velocity. Thus, if the mass of the dust disk 
is dominated by 10 /zm particles (for example, if 
the particles have size distribution n oc s~^'^ with 
maximum size 10 /^m), the dust-gas ratio would 



increase through the gas accretion. For example, 
if the initial mass of the gas disk is O.IIMq and 
it reduces to O.OIM© after viscous accretion, the 
dust-gas ratio would increase to be 6 times of its 
initial value. The increase in the dust-gas ratio 
speeds up the formation of planetesimals through 
mutual collisions. Their enhanced abundance may 
also cause the dust component to become unstable 
and promote the planetesimal formation through 
the gravitational instability (Goldreich & Ward 
1973; Sekiya 1998). 

If the particle growth proceeds to make 1 mm 
to 1 cm particles during the gas accretion phase, 
such large particles migrate inward rapidly and 
accumulate in the inner part of the disk (see Fig 
9a). This size fractionation causes an increase in 
the dust-gas ratio at the inner disk, while this pro- 
cess may decrease the dust-gas ratio at the outer 
disk, resulting in a dust disk concentrated to the 
inner part of the gas disk. 

4.3. Summary 

The radial migration of dust particles in accre- 
tion disks is studied. Our results arc as follows. 

1. Dust particles move radially both inward and 
outward by the gas drag force. Particles at high 
altitude (|z[ > 2hg) move outward because they 
rotate slower than the gas whose pressure gradient 
force is inward. Small particles (s < 100 fim at 
10 AU) near the midplane (|^;| < hg) are advected 
by the gas outflow. On the other hand, particles 
at intermediate altitude and large particles (s > 
1 mm at 10 AU) move inward. 

2. The net radial velocity, averaged in the verti- 
cal direction, is usually inward, provided that the 
radial gradient of the gas surface density is not too 
steep {ps > —1.3). Sedimentation removes out- 
flowing particles from high altitudes. Small par- 
ticles, which can be advected by the outflowing 
gas around the midplane, do not concentrate at 
the midplane. In the inner part of the gas disk 
(r < 100 AU for 10 /im particles), the inflow ve- 
locity of particles is smaller than the gas accretion 
velocity, resulting in an increase in the dust-gas 
ratio. 

3. The particle sedimentation would be efficient 
if dust-gas coupling is relatively weak (Sc > 1.0). 

If the sedimentation is so efficient (Sc > 10), the 
number of outflowing particles around the mid- 
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plane would be large, and the direction of the net 
radial velocity of particles would change to out- 
ward at some distances from the star. Accumula- 
tion of particles at such locations serves to increase 
the local dust-gas ratio. 

4. The inflow velocity of particles depends on 
the particle size. Therefore, the inflow causes the 
size fractionation of particles. Larger particles ac- 
cumulate at distances closer to the star. 
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Fig. 1. — Radial velocity of gas Vr^g. The distance 
from the midplane z is normalized by the disk scale 
height hg. In the standard model {q = —0.5), the 
radial velocity is independent of the distance from 
the star. The dotted line shows the zero velocity 
for reference. 
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Fig. 3. — Radial velocity Vr,d of dust particles of 
s = 10 /Ltm (solid line), 100 /im (dashed line), and 
1 mm (dot-dashed line) at 10 AU. The distance 
from the midplane z is normalized by the disk scale 
height hg = 0.59 AU. The dotted line shows the 
radial velocity of the gas Vr^g. 



Fig. 2. — The particle radius whose non- 
dimensional stopping time Tg is unity in the stan- 
dard disk {ps = — 1.0 and q = —0.5). The solid line 
corresponds to particles at the midplane, while 
the dashed line corresponds to particles at altitude 
z = 2hg. Particles smaller than the lines have Tg 
smaller than unity and couple well with the gas, 
while particles larger than the lines have > 1. 
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Fig. 4. — Dust-gas ratio, PdJ pg-, against the dis- 
tance from the midplanc (o) for various sizes at 10 
AU and (6) at various locations for 100 /^m. The 
distance from the midplane z is normalized by the 
disk scale height hg. The averaged dust-gas ratio 
is assumed to be 0.01. (a) Three solid lines show 
the dust-gas ratio of particles of 1 mm, 100 /zm, 
and 10 /im from the left lino. (6) Three solid lines 
show the dust-gas ratio at 100 AU, 10 AU, and 1 
AU from the left line. 
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Fig. 5. — Net radial velocity of particles well 
mixed with the gas, {vr,d)- The three solid lines 
correspond to s = 10 /im, 1 /im, and 0.1 /im par- 
ticles from the upper line. The dashed line shows 
the gas accretion velocity. 
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Fig. 6. — Normalized mass flux pdVr,d/^d at 10 
AU. The solid, dashed, and dot-dashed lines are 
for s = 10 /Ltm, 100 /xm, and 1 mm particles, re- 
spectively. For 1 mm particles (dot-dashed line), 
1/10 of the value is plotted (or refer the ordinate 
at the right hand side). The distance from the 
midplane z is normalized by the disk scale height 
hg = 0.59 AU. The dotted line shows the zero ve- 
locity for reference. 
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Fig. 7. Net radial vcloc;ity {vr,d) of scdimcntcd 
particles. Three solid lines correspond to s = 1 
mm, 100 yum, and 10 /xm particles from the left 
line. The accretion velocity of the gas is shown 
by the dashed line, and the dotted line shows the 
zero velocity for reference, (o) the standard model 
(p^ = -1.0, q = -0.5, and Sc = 1.0) {b) = 
-0.5. (c) = -1.3. (d) Sc = 10. 
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Fig. 8. — Function F for various values of the 
Schmidt number, Sc = 10, 1, and 0.1 from the 
upper line. The dotted line shows the zero for 
reference. 
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Fig. 9. Surface density of dust disks, S^. (a) 
The disk is composed of single sized particles of 
s = 1 mm, 100 /Ltm, and 10 /xm from the left line. 
The mass flux M^i is assumed to be 10~"'^''Mq yr"-'^. 
(The mass flux of the gas is 10"^Mq yr'^) The 
gas disk is the standard model {ps = —1.0). (6) 
The disk is composed of particles with a power 
law size distribution. The indices of the surface 
density profiles of the gas disks are Ps = —1.0 
and —0.5 from the upper line. The particle num- 
ber flux is assumed to be proportional to s~^-^. 
The minimum and maximum sizes of particles are 
0.1 /im and 10 cm, respectively. The total mass 
flux is Md.aii = 10-i°Mq yr-i. 
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